Energy loss and thermalization of heavy quarks in a 
strongly-coupled plasma 



C. Marquet", G. BeuP and B.-W. Xiao* 

" lnstitut de Physique Theorique, CEA/Saclay, 91191 Gif-sur-Yvette cedex, France 
b Nuclear Science Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 



Abstract 

Using the AdS/CFT correspondence, we compute the medium-induced energy loss of a de- 
celerating heavy quark moving through a strongly-coupled supersymmetric Yang Mills plasma. 
In the regime where the deceleration is small, a perturbative calculation is possible and we ob- 
tain the first two corrections to the energy-loss rate of a heavy quark with constant velocity. The 
thermalization of the heavy quark is also discussed. 



1. Motivation 

Nuclear modification factors Raa compare particle production in A+A collisions to a set 
of independent p+p collisions, such that Raa is unity in the absence of medium effects. In 
relativistic heavy ion collisions, in which a thermalized quark-gluon plasma (QGP) is formed, 
the suppressed production of high- pr particles (Raa < 1) is interpreted as a consequence of the 
interaction of the plasma with the hard partons that initiate the hadrons. Indeed, after an energetic 
parton is produced in the early stages of the collision, it loses energy while propagating through 
the QGP, and may even thermalize before exiting the medium, meaning less particle production 
at high-p r . 

When one tries to reproduce the value Raa — 0.2 observed for light hadrons with perturba- 
tive QCD (pQCD) medium-induced energy loss calculations, one obtains that the jet quenching 
parameter q, which essentially characterizes how fast energetic partons lose energy, is of order 
10 GeV 2 /fm fH, a value unnaturally high for a weakly-coupled QGP. Moreover, once this pa- 
rameter is fixed, pQCD calculations predict that the Raa for heavy hadrons should be higher than 
it is for light hadrons, since heavy quarks lose energy at a smaller rate than light quarks. The fact 
that the Raa $ for light and heavy hadrons are comparable adds to the puzzle and one wonders 
whether the pQCD framework is applicable in this context. 

This motivated to think about strongly-coupled plasmas. Addressing the strong-coupling 
dynamics in QCD is an outstanding problem, lattice simulations remain the only method to 
obtain quantitative results, but are inefficient when real-time dynamics is required, which is the 
case when analyzing jet quenching. However for a class of non-abelian thermal gauge theories, 
the AdS/CFT correspondence |3, l2l |4J] provides an alternative, and it has become a popular 
approach to assume that just above the critical temperature T c where the conformal anomaly in 
QCD is small, the QGP can be approximated by an N — 4 Super Yang Mills (SYM) plasma in 
the strong coupling regime. Information about real-time dynamics can be obtained within this 
setup, and in this work we address the energy loss and thermalization problem for heavy quarks 
propagating through the SYM plasma. 
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2. The trailing-string picture 



We consider the large N c , small gauge coupling gy/u limit of the SYM theory: 

N c -> oo , gYM -> , /I = g\ M N c finite, (1) 

where the 't Hooft coupling A controls the theory. Then strong coupling means A » 1. In this 
regime the equivalent string theory in AdS 5 space is weakly coupled and weakly curved: 

gYM < 1 <=> g s «: 1 and A » 1 /? » Z s , (2) 

where g, is the string coupling, l s is the string length and R is the curvature radius of the AdS 
space. In this limit of small string coupling and large curvature radius, classical gravity is a good 
approximation of the string theory, and the background metric corresponding to the thermal SYM 
theory is 

r du 2 

ds 2 =R 2 h(u)dt + u 2 (dx') 2 = GuydX^dX* (3) 

h(u) 

with h(u) — u 2 - ui/u 2 and where m/, = nT is a black-hole horizon. The corresponding Hawking 
temperature T is the temperature of the SYM plasma. The coordinate in the fifth dimension 
m = r/R 2 has the dimension of momentum and the SYM theory lives on the boundary at u = 00. 

The heavy quark whose energy loss we want to compute is in the fundamental representation, 
and lives on a brane at u — u m — 2mn/ yfl — > 00 with m the vacuum mass of the quark. A string 
is attached to it and hangs down towards the horizon. Points on the string can be identified to 
quantum fluctuations in the heavy quark wave function with virtuality ~ u. Indeed, the quantum 
dynamics in the SYM theory is mapped onto classical dynamics in the 5th dimension. More 
precisely, the string dynamics is given by the Nambu-Goto action 



2nR 2 



j ' dTda-^Pg with g = detg a 



where r and cr are the worldsheet coordinates and g a t, = Gfj V (d a X M )(dhX v ) is the metric induced 
on the worldsheet. Let us assume that the quark moves along the x direction and parameterize 
the space-time coordinates in the following way: X^ = (t, x, y, z, u) = (r, x(r, cr), 0, 0, cr). Then 

yfA r I u 2 * 2 " ~ 

S — I dt du -ill h u 2 h(u)x' 2 (5) 

In J ^ h(u) 



with x = d t x(t, u) and x' - d u x(t, u). The classical equation of motion d a 5£,/5d a X l = gives 
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From the solution of this equation, one gets the rate at which the energy flows down the string: 

C"-£ yA. R 2 

m) = 7F-^o = - =u 2 h{u)xx! . (7) 

In JUS], the authors imagined using an external force to compensate the energy loss and pull 
the quark at a constant velocity v. Writing x(t, u) — xq + vt + vF(u), they obtained 

(8) 
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In this stationary setup, the corresponding energy-loss rate (0 is also independent of u and is 
therefore identified to the heavy-quark energy loss: 



dE ^22 fQ , 

Investigating this picture in more detail, one finds that the energy loss is due to the presence of an 
induced horizon on the string worldsheet [7]. This special point in the fifth dimension, located at 
u — u s where 

u s = Vr«/! with y = , (10) 

divides the trailing string into two parts. This is the location where the local speed of light in AdS 
space coincides with the propagation speed of the string: no light ray emitted from the part of the 
string below u s can reach the heavy quark (even if it will reach the boundary at u m ), the part of the 
string with u < u s is not causally connected to the quark. In |8], it was shown that the worldsheet 
horizon corresponds to the saturation scale Q s on the gauge theory side (the scale at which the 
scattering of dilute probes off the plasma becomes strong), and that the quantum fluctuations 
corresponding to points on the string with u < u s should be thought of emitted radiation, rather 
than part of the heavy-quark wave function. The picture is identical for a weakly-coupled plasma, 
but with a different saturation scale |9|]. 

3. Letting the heavy quark decelerate 

In the following we study the system formed by the upper part of the string, dual to the heavy 
quark dressed with quantum fluctuations. The energy of this system is 

f "»< 6£ yfAR 2 r u " 1 + u 2 h(u)x' 2 

Using the equations of motion one can derive the following energy-conservation equation ifioll : 

dE 6 E\ du 

— = — — - - <D(f, u s ) + <D(f, u m ) , (12) 
at at 

where 5 E\ du jdt is the energy change of the system due to possible variations of u s . The energy 
flow <t>(f , Um) represents the work external forces acting on the heavy quark while the energy flow 
<l)(f, u s ) represents the energy radiated into the plasma by the dressed quark. 

In 01, an external force is acting on the heavy quark and adjusted in order to have a stationary 
solution with dE/dt = 6 E\j u /dt - 0. The energy flow is also uniform along the string, agreeing 
with 0(f, u s ) = Um) : the same amount of energy is put into the string at the top than flows 
at the bottom into the worldsheet horizon, or equivalently is radiatively lost by the heavy quark. 
In this work we would like to consider the case where the external force is turned off and the 
heavy quark is allowed to decelerate. Indeed using <]>(?, u m ) = in (fT2l imposes that — v oc u, ju m . 
Assuming a slow deceleration, the energy loss for this setup can be computed perturbatively. 
Although the total energy loss -dE/dt of the upper part of the string includes the contribution 
from the term 5 E\ dUs , we focus on the radiative energy loss u s ) which is physically more 
important as it is directly related to the deceleration of the quark, and therefore control how fast 
thermalization is achieved. 
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With first-order (v) and second-order (v 2 and v) corrections, inserting the Ansatz 
x(t,u) = x + j v(Orff' + v(f)F( M ) + v(«) ((k, v(f)) + i' 2 (0,ri («, HO) + v(0 ^2 v(0) (13) 



into the equation of motion (O, one obtains the first-order deformation of the trailing string: 
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This allows to compute the energy flow <J>(f, u), and imposing 0(f, u m ) = yields 



m = - 



y 2 (u m - u s ) 
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(14) 



(15) 



and one can now specify our expansion parameter: m/,/(m ot - u s ). In particular our perturbative 
expansion breaks down if u s is too close to u m , meaning if the worldsheet horizon is too close to 
the flavor brane. Finally the radiative energy loss is 



<£(f, u s ) = ^-yv 2 u 2 h 
Zn 



1 - 



(1 + y 2 )u 2 F(u s ) 
y 2 (u m - u s ) 
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(ll m Us)^ 



(16) 



As F (u) is negative, the first-order correction is positive, meaning that the rate of energy loss gets 
stronger with time. However this is only efficient for small v because in the ultrarelativistic case 
y » 1, the first-order correction is suppressed by a factor y^l 2 . As v decreases, the string slowly 
gets straighter to eventually become that of a static quark but we do not expect our picture to 
apply that far. Before this happens the heavy quark will thermalize and the discussion should then 
be modified to include the transverse Brownian motion of the quark. However our calculation 
allows to estimate how fast thermalization is reached. 

Expressions including second-order corrections can be found in [101. For the energy loss, this 
contribution corresponds to radiation due to the quark deceleration. For dv/dt, one can write: 
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(17) 



(p " Vt) 2 

where t — nT t and yu = m,„/m/, = 2m/ ( yfAT) (note that /i » ^Jy in order for our perturbative 
expansion to be valid). We recover the estimate of the thermalization time made in ylla]: 

r=^~, d8) 
but our corrections lead to thermalization faster than what is inferred from the stationary solution. 
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